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^ 1 ABSTRACT 

O i, We study the impact of the environment on the evolution of galaxies in the zCOSMOS 10k sample in the redshift range 0.1 < z < 1.0 over an area 
I of ~ 1.5 deg 2 . The considered sample of secure spectroscopic redshifts contains about 8500 galaxies, with their stellar masses estimated through 
O ' SED fitting of the multi-wavelength optical to Near-InfraRed photometry. The evolution of the Galaxy Stellar Mass Function (GSMF) in high and 
,£j , low density regions provides a tool to study the mass assembly evolution in different environments; moreover, the contributions to the GSMF from 
t/^ . different galaxy types, as defined by their SEDs and their morphologies can be quantified. At redshift z ~ 1 the GSMF is only slightly dependent 
, ^ , ■ on environment, but at lower redshifts the shapes of the GSMFs in high- and low-density environments become extremely different, with high 
density regions exhibiting a marked bimodality, not reproducible by a single Schechter function. Galaxy evolution depends both on the stellar mass 
I , and on the environment, setting the probability of a galaxy to have a given mass. All the galaxy properties related to the stellar mass will show a 
, dependence on the environment, reflecting the difference observed in the mass functions. The shape of the GSMFs of early and late type galaxies 
• is quasi-universal, i.e., it is consistent to be the same in the extreme environments we considered. The evolution toward z = of the transition 
I ' mass Alcross, i.e., the mass at which the early and late type GSMFs match each other, is more rapid in high density environments. The same 
result can be obtained studying the relative contribution of different galaxy types, implying a faster evolution in overdense regions, in particular 
f~>) , for intermediate stellar masses. The rate of evolution is different when dividing galaxy types on the basis of their SEDs or their morphologies, 
• 1 tentatively suggesting that the migration from the blue cloud to the red sequence is occurring on a shorter timescale than the transformation from 
C ' disc-like morphologies to ellipticals. Environmental mechanisms of galaxy transformation are at work at z < 1. The comparison of the observed 
| GSMFs to the same quantities derived from a set of mock catalogues based on semi-analytical models shows a disagreement, both in low and high 
ON . density environments. 
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$H ' 1 . Introduction sity environments coincided with galaxy clusters or groups and 

the more general effect of the environment on field galaxy evolu- 

The environmental dependence of galaxy properties (colour, star tion is still poorly explored . The evolution of the Galaxy Stellar 

formation, mass) is a well known effect in the local universe. At Mass Function (GSMF) as a function of the large scale environ- 

present many local studies have been carried out to analyse the ment has bee n studied in the DEEP2 Galaxy Redshift Survey 

influence of environment on colours, luminosities, morpholo- (Bun dv et al.ll200rjh . considering the redshift range z = 0.4 - 1 .4, 

gies, structural parameters, star formation and stellar masses: all which Umits me connec tion between this study and those in the 

local relations are somewhat link ed to the morphology-density i oca i Universe 

relation shown by Dressier ( 1980). 

. i . . , ... c . , , Some of the still open question are: which is the leading 

At higher redshifts this kind of study becomes very ditti- .. * nT * ^ , 

, ■ c , • iff-, property in the evolution of field galaxies Is the fate of a galaxy 

cult, because the need of large spectroscopic samples of taint . ■ , ^ , , , , i 

. . ... , i f. i- li decided once its mass is defined or do some external players have 

galaxies with a good sampling rate hampers a reliable estimate , n . , ... . . , , , , , 

j.,, y t i ~ r c , , a role And, it the environment plays such a role, when does it 

of the environment. Up to now, therefore, most of the high den- „ , , . j , , ■ • , , • n 

start to affect galaxy evolution, and through which mechanism ! 
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e-mail: micol .bolzonellaOoabo . inaf . it galaxies is not unanimous. At low redshifts. lBalogh et alJd2001l) , 

* Based on observations obtained at the European Southern usin 8 2MASS and LCRS data, separated different environments 

Observatory(ESO) Very Large Telescope (VLT), Paranal, Chile, as as field, groups and clusters, finding that Luminosity and Mass 

part of the Large Program 175.A-0839 (the zCOSMOS Spectroscopic Functions depend both on galaxy type (with steeper functions 

Redshift Survey). for emission line galaxies) and environment (with more massive 
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and brighter objects more common in clusters), mainly as a con- 
sequence of the different contribution of passive galaxies. 

Most of the low-redshift studi es are based on the SDSS data: 
studying the colour bimodality, Balog rTet al.l (|2004) found that, 
even if the mean colour of red and blue galaxies depends only 
marginally on environment, the fraction of red galaxies signifi- 
cantly increases toward high densities at fixed luminosity; they 
propose that the properties of star forming galaxies are mainly 
related to their mass and that, to preserve the bimodality with- 
out altering the colours modelled by two Gaussian distributions, 
the transformation from late to early type galaxies should be 
rapid in truncating the star formation and efficient for all lu- 
minos ities and environments. In an analogous study. [Hogg et al.l 

(2003) reached similar conclusions, with blue galaxies showing 
no correlation between their luminosity/mass and local density 
at a fixed colour. In a comprehensive analysis o n environmental 
depen dence of galaxy properties in the SDSS, Kauffmann et alj 

(2004) found evidences supporting a mix of both nature and nur- 
ture scenarios, where the star formation history is the galaxy 
property more dependent on environment. Tanakaet al. (2004) 
studied the environmental dependence of star formation and 
morphologies: they pointed out that environment is more effec- 
tive on faint galaxies, affecting more their star formation and 
morphology in dense regions above a critical density, whereas 
the properties of bright gala xies do not show any dependence on 
environment. iBlanton et al.l ([2005 ) found that the star formation 
history and the mass of bright galaxies of the SDSS are more 
dependent on environment than struct ural paramet e rs, lik e the 
surface brightness and profile shape. Bal drv et"al1 (120061) . us- 
ing a sample limited at z < 0.1 and a projected density esti- 
mator, found that the characteristic mass contributing the most 
to the stellar mass density increases with density from 10 10 6 to 
10 10 9 A1 Q and the fraction of galaxies in the red sequence in- 
creases as a function of mass and density: they concluded that 
both the stellar mass and the environment, which cannot be con- 
sidered as a second order effect, determine the prob ability of a 
galaxy being in the red sequence. Bal drv et"aTI ((2008) computed 
the local estimate of the GSMF detecting a low mass upturn at 
~ 10 9 A4 Q , like the one observed in galaxy clusters, but the dom- 
inant galaxy types and the processes producing the bimodality 
of the GSMFs shou ld be different in clusters and in the field. 
Mate us et al.l {2008) used SDSS data and compared their results 
to the Millennium simulation: they also supported the nurture- 
like scenario, with environment playing a fundamental role in 
defining the galaxy properties and confirmed the findings by 
Mate us et al.l J2007): galaxy evolution takes place in an acceler- 
ated way at high redshift in denser environments, wher e massive 
galaxi es preferentially reside. Following the study by iLee & Lil 
(2008) on SDSS, the initial conditions significantly contribute 
in establishing the environmental dependence of galaxy proper- 
ties, even if in the densest regions the initial correlation between 
galaxy properties and environment will be superseded by non 
linear processes at small-scales. Recently. Ivan den Bosch et al] 
(2008) investigated the efficiency of transformation processes, 
using a group catalogue drawn from the SDSS: the colour and 
the concentration of a satellite galaxy are almost completely de- 
termined by their stellar mass, with only a very mild dependence 
on environment. 

At higher redshifts, probing the effect of environment on 
galaxy evolution becomes more difficult and often this kind of 
studies uses projected estimators of local density and relies on 
photometric redshifts. The main spectroscopic study carried out 
up to now estimati ng the effect o f envi ronment on GSMFs has 
been presented by Bundy et al. (2006): using DEEP2 data at 



0.4 < z < 1 .4 and Rab < 24. 1 they drew the conclusion that the 
quenching of star formation, and then the transition between the 
blue cloud and the red sequence, is primarily internally driven 
and dependent on mass, even if they detected a moderate ac- 
celeration of the downsizing phenomenon in overdense regions, 
where the rise of the quiescent population with cosmic time ap- 
pears to be faster, as seen through the evolution of the transition 
and quenching masses, A1 croS s and A1q. 

Us ing the same dataset complemented by SDSS at low red- 
shifts, ICooper et alj d2008l) studied the connection between the 
Star Formation Rate (SFR) and environment, finding hints for 
a reversal of that relation from z ~ 0, where the mean SFR 
is decreasing with local density, to z ~ 1, where a blue pop- 
ulation causes an increase of the mean SFR in overdense re- 
gions; nonetheless, the decline of the global cosmic star for- 
mation history since z ~ 1 seems to be due to a gradual gas 
consumption rather than environment-dependent processes: the 
fall of the mean SFH from z ~ 1 to z ~ is much stronger 
than its dependence on environment, which can be regarded as 
a small perturbing effect to the global trend. A similar result on 
th e reversing relatio nship SFR-environment has been achieved 
bv lElbaz et al]d2007l) , using GOODS data and SFR derived from 
UV and 24 //m e mission. Looki n g at th e evolution of the colour- 
density relation, ICucciati et al.l (|2()06) used spectroscopic data 
from the VVDS up to z ~ 1.5: they found a steep relation at 
low-z, tending to disappear at higher redshifts. In particular, they 
identified differences in colour distributions in low and high den- 
sity regimes at low redshifts, whereas at high redshifts the envi- 
ronment is not affecting such a distribution. In their proposed 
scenario the processes of star formation and gas exhaustion are 
accelerated for more luminous objects and high density environ- 
ments, leading to a shift with cosmic time of the star formation 
activity toward f ainter galaxies and low density environments. 
IScodeggio et al.l ((2009) studied the stellar mass and colour seg- 
regations in the VVDS at redshifts z = 0.2 - 1 .4, using a density 
field computed on scales of ~ 8 Mpc; they found that the colour- 
density relation is a mirror of the stellar mass segregation, that in 
turn is a consequence of the dark matter halo mass segregation 
predicted by hierarchical models. 

Both the effects of environment on local galaxy properties 
and on their evolution are still uncertain, making the nature vs 
nurture debate still open. From the aforementioned results, there 
seems to be some hint that the galaxy evolutionary path from 
the blue cloud to the red sequence depends on environment, 
but the determination of the mechanism of such a transforma- 
tion, its probability, its link to the environment and to the innate 
galaxy properties is a difficult task. Different physical processes 
of galaxy transformation will differ on timescales, on efficiency 
and on the observational repercussions, like colour and morphol- 
ogy. The GSMF is a very suitable tool to investigate this problem 
and to witness the buildup of galaxies and its dependence on en- 
vironment. 

In this paper we will focus on the effect of environment on 
field galaxies, avoiding the most extreme overdense regions like 
cluster cores, and we will use data from COSMOS (Cosmic 
Evolution Survey) and zCOSM OS; parallel and complemen- 
tary analyses are pre s ented in Pozzetti et al. (2009), Zucca et al. 
(120091) , llovino et all d2009l) . ICucciati et all (120091). iTasca et al. 
d2009l) . Kovac et al. d2009l) and Vergani et al. d2009l) . The plan 
of this paper is the following: in Sec. [2] we describe the spectro- 
scopic and photometric datasets and the derived properties we 
used to characterise different galaxy populations; in Sec. [3] we 
derive the GSMFs and in Sec. [4] we analyse the different con- 
tribution of galaxy types to the GSMF in different environment, 
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discussing the implications for the picture of galaxy evolution in 
Sec.0 

Throughout the paper we adopted the cosmological parame- 
ters Q m = 0.25, Q A = 0.75, H = 70kms _1 Mpc 1 and magni- 
tudes are given in the AB system. 

2. Data 



The zCOSMOS survey dLillv et al.ll2007l) is a redshift survey in- 
tended to measu re the distances of ga laxies and AGNs over the 
COSMOS field dScoville et al.ll2007bh . The whole field of about 
2 deg 2 has been observed from radio to X-ray wavelengths by 
parallel projects, involving worldwide teams and observatories. 
The starting point has been the largest HST survey carried out up 
to now with ACS. The coexistence of multi-wavelength obser- 
vation, morphologies, spectroscopic redshifts make COSMOS 
a unique opportunity to study the evolution of galaxies in their 
large scale structure context. 

2.1. Spectroscopy 

The spectroscopic survey zCOSMOS is currently on-going and 
is subdivided in 2 different parts: the "bright" survey, aiming 
at observing ~ 20000 galaxies, with a pure flux limited selec- 
tion corresponding to 15 < TAB < 22.5, and the "deep" survey, 
whose goal is the measure of redshifts in the range 1 .4 < z < 3.0, 
within the central 1 deg 2 . 

The data used in this paper belong to the so-called 10k sam- 
ple dLillv et al.ll2009l) . consisting of the first 10644 observed ob- 
jects of the "bright" survey, over an area of 1.402 deg 2 with a 
mean sampling rate of ~ 33%. The final design of the survey 
aims at reaching a sampling rate of ~ 60 - 70%, achieved by 
means of an eight-pass strategy. The observations have been 
carried out with VIMOS@VLT with the red grism at medium 
resolution R ~ 600. T he data have been reduced with VIPGI 
dScodeggio et alj 120051) and spectroscopic redshifts have been 
visua lly determined after a first hint provided by EZ dGarillil 
20090 The reliability of the redshift measurements has been 
represented by means of a flag: a 4 value corresponds to 100% 
confidence, 3 to ~ 90%, 2 to ~ 75%, 1 to ~ 50% and to a non 
measured redshift. All the redshifts have been checked by more 
than two astronomers. A decimal digit spec ifies if the redshift is 
in agreement with photometric redshifts (Feld mann et al.i r2006) 
computed from optical and near infrared (NIR) photometry us- 
ing the code ZEB RA (Zurich E xtragalactic Bayesian Redshift 
Analyzer, Feldmann et al. 2008). Moreover for some objects the 
measure resulted to be hampered by technical reasons (for in- 
stance the spectrum at the edge of the slit); in those cases a 
flag -99 has been assigned. Different flags have been assigned 
to identify broad line AGNs and targets observed by chance in 
slits. 

2.2. Photometry 

The photometry used in the following is part of the COSMOS 
observations and encompasses optical to NIR wavelengths: u* 
and K s from CFHT, Bj, Vj, g + , r + , i + and z + from Subaru, and 
Spitzer IRAC magnitudes at 3.6, 4.5, 5.8 fim. Detai ls on photo- 
metric observations and da t a redu ction are given in Capa k et alj 
d2007bl) : iMcCracken et alj d2009l) . The scantiness of standard 
stars in photometric observations and the uncertainty in the 
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knowledge of the filter responses result in an uncertain calibra- 
tion of zero-points. To avoid this inconvenience, we optimized 
the photometry by applying offsets to the observed magnitudes: 
we computed such photometric shifts in each band minimizing 
the differences between observed magnitudes and reference ones 
computed from a set of Spectral Energy Dis tributions (hereafter 
SEDs). We adopted an approach similar to Cap ak et all ([2007b, 
see their Table 13), but considering the same set of SEDs we 
used to compute stellar masses detailed in Sec. 12.31 obtaining in 
general very similar offsets for all the filters. 

2.3. Stellar masses 

Stellar masses have been evaluated by means of a SED fitting 
technique, using the code Hyperzm ass, a modified version of the 
photometric re dshift code Hyperz ([Bolzonella et aU uOOO). In a 
recent paper, Marchesini et al.1 d2008l) analysed the effect of ran- 
dom and systematic uncertainties of the stellar mass estimates on 
the GSMF, considering the influence of metallicity, extinction 
law, stellar population synthesis model, and initial mass func- 
tion. Here we describe the approach and the tests we performed 
on our data. 

We used different libraries of SEDs, derived from differ- 
ent models of stel l ar pop ulation synthesis: 1. the w ell-known 
iBruzual & Charlotl d2003l here after BC03) library , 2. iMarastonl 
d2005l hereafter M05) and 3. ICharlot & Bruzuall d2009j here- 
after CB07). The main difference among the three libraries is 
the treatment of thermally pulsing asymptotic giant branch (TP- 
AGB) stars. M05 models include the TP-AGB phase, calibrated 
with local stellar populations. This stellar phase is the dominant 
source of bolometric and near-IR energy for a simple stellar pop- 
ulation in the age range 0.2 to 2 Gyr. Summing up the effects of 
both overshooting and TP-AGB, the M05 models are brighter 
and redd er than the BC03 mo dels for ages between ~ 0.2 and 
~ 2 Gyr dMaraston et alj|2006l) . The use of the M05 models im- 
plies the derivation of lower ages and stellar masses for galaxies 
in which the TP-AGB stars are significantly contributing to the 
observed SED (i.e., ages of the order of ~ 1 Gyr). At older ages 
the M05 models are instead bluer. CB07 is the first release of 
the new version of the Charlot& Bruzual li brary, not yet pub- 
lic. CB 07 models include the prescription of iMarigo & Girardil 
(2007) for the TP-AGB evolution of low and intermediate-mass 
stars. Similarly to M05 models, this assumption produces signif- 
icantly redder near-IR colors, and hence younger ages and lower 
masses, for young and intermediate-age stellar populations. A 
brief description of the effect on GSMFs induced by different 
ch oices of template SED s can be found in the companion paper 
bv lPozzetti etafl d2009l) . 

All the considered libraries provide a Simple Stellar 
Population (SSP) and its evolution in many age steps for a fixed 
metallicity and a given Initial Mass Function (IMF); it is possible 
from the SSP models to derive the composite stellar populations 
that can reproduce the different types of observed galaxies, im- 
posing a star formation history (SFH). We built 10 exponentially 
declining SFHs with e-folding times ranging from 0. 1 to 30 Gyr 
plus a model with constant star formation. Smooth SFHs are a 
simplistic representation of th e complex SFHs gal axies can have 
experienced. Nonetheless in Pozzetti et al. (2007) we also com- 
puted stellar masses using SEDs with random secondary bursts 
superimposed to smooth SFHs, finding differences well within 
the statistical uncertainties for most of the sample. Instead, for a 
minor subsample of galaxies, stellar masses obtained using mod- 
els with secondary bursts can become ~ 40% larger, although the 
effect on the final GSMF remains negligible. 
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The IMF is another important parameter: differences be- 
tween stellar mass estimates derived with different underly- 
ing IMFs can be s tatistically re covered. T he most used IMFs 
are the S abeter dSalpeter||1955l) . Kroupa dKroupal |200T1) . and 
Chabrier (Chabrier 2003) ones. The statistical differences of 
stellar masses is logM^ ~ lo S^ciiab + °- 23 ' ^S^Chab ~ 
log Mt - 0.04. Besides the mass, the best fit SEDs obtained 
with 2 different IMFs can correspond to significantly different 
values for other parameters used in the SED fitting process. 

In stellar population synthesis models the metallicity can be 
evolving or fixed. We have verified in simulated catalogues that 
introducing different values of metallicities does not introduce a 
substantial improvement of the quality of the best fits, at the cost 
of the introduction of an additional parameter. Consequently we 
adopted a fixed solar metallicity, even though this choice may 
introduce another small uncertainty on the stellar mass estimate. 

Dust extinction has b een modeled by means of the Calzetti's 
law dCalzetti et al.ll2.Q00b . with values ranging from to 3 mag- 
nitudes of extinction in V band. 

The x 1 minimization comparing observed and template 
fluxes at a fixed redshift z - z S pec provides the best fit SED, to 
which are associated a number of physical parameters, like the 
age, the reddening, the instantaneous star formation and the stel- 
lar mass. Notice that the meaning of stellar mass throughout this 
paper is not the integral of the star formation, because from that 
value we have to exclude the return fraction, i.e., the fraction of 
gas processed by stars and returned to the interstellar medium 
(ISM) during their evolution. 

Tests on simulated catalogues considering the effect on stel- 
lar mass estimates of different choices of reddening law, SFHs, 
metallicities, SED libraries show a typical dispersion of the or- 
der of ( TioeM - 0.20. Ev en a simpler technique lik e the one 
used b y iMaier et aTl (120091) and derived from Eq. 1 of iLin et all 
(2007j), produces a scatter not larger than ~ 0.16 dex, although 
with some small trend as a function of stellar mass and redshift. 
These tests show that stellar mass is a rather stable parameter 
in SED fitting when dealing with a set of data spanning a wide 
wavelength range extended up to NIR. 

Since the fluxes provided by the available libraries at IR 
wavelengths are extrapolated, the choice of filters used in best 
fitting is reduced to 2.5 yum rest-frame for M05 models (at longer 
wavelengths these models are using the Rayleigh-Jeans tail ex- 
trapolation) and to 5 fim rest-frame for BC03 and CB07 models, 
since at longer wavelengths also the dust re-emission can con- 
tribute to the flux budget. 

A problem arising when dealing with a very large number 
of template SEDs is to av oid non-physical bes t fits. We applied 
2 priors (the same u sed inJPozzetti et al. 2007, and proposed by 
iFontana et al.ll2004l and fKauffma nn et al.ll2003l) to avoid such a 
problem. In particular, we excluded best fit SEDs not fulfilling 
the following requirements: 1. Ay < 0.6 if age/r > 4 (i.e., old 
galaxies must have a moderate dust extinction); 2. star formation 
must start at z > 1 if r < 0.6 Gyr. Moreover, we tested through 
simulations that imposing a minimum best fit age of 0.09 Gyr 
reduces potential degeneracies and improves the reliability of 
the stellar mass estimate. The maximum allowed age is the age 
of the Universe at z spe c. 

As mentioned in Sec. 12.21 the first SED fitting run over the 
brightest galaxies and most secure galaxy redshifts has been per- 
formed to compute the photometric offsets. We checked that fur- 
ther iterations of the SED fitting and offset estimation do not 
improve significantly the^- 2 statistics. 

To ease the comparison with literature results, in the fol- 
lowing we present GSMFs obtained adopting the BC03 stellar 



masses. However the qualitative trends are the same for any 
choice of the stellar population synthesis models. 

2.4. Environment 

The density field has been derived for the 10k spectroscopic 
sample using different estimators co mbined with the ZADE 
(Zurich Adaptive Density Estimator, iKovac et"a"D |2009) algo- 
rithm. Most of the existing studies rely on photometric redshifts 
and projected densities computed in wide redshift slices, pos- 
sibly diluting the signal from overdense regions. An important 
added value of COSMOS is the availability of spectroscopic red- 
shifts obtained with a good sampling rate, making feasible an 
accurate estimate of the environment, with high resolution also 
on the radial direction. 

To this aim, spectroscopic redshifts have been used as a 
skeleton of galaxy structures, with the incorporation of a sta- 
tistical treatment of the likelihood function of photometric red- 
shifts computed with ZEBRA. This approach allows to probe 
a wide range of environments and to reduce the Poisson noise. 
Results have been extensively and carefull y tested on mock cat- 
alogu es from the Millennium simulation (Kitz bichler & Whitel 
2007). The reconstruction of overdensities 1+6 has been ex- 
plored using different tracer galaxies, different spatial filters and 
different weights (e.g., luminosity or stellar mass) assigned to 
each galaxy. The density contrast 5 is defined as (p-p)/p, where 
p is the density as a function of RA, DEC and z and p is the 
mean density measured at the same redshift. In principle, the 
fully realistic physical representation of the environment should 
involve the mass of the dark matter haloes in which the galaxies 
are embedded. This mass is clearly not directly accessible to ob- 
servations, hence an affordable surrogate to weight the number 
density field is given by the stellar masses of the surrounding 
galaxies. This is a proxy of the overall density field, since galax- 
ies are biased tracers of the underlying matter distribution, and 
stellar mass lacks the gas component that differentiates red and 
blue galaxies. Since the choice of a fixed selection band results 
in different populations preferentially sampled at different red- 
shifts, weighting with stellar mass should also mitigate this is- 
sue. As expected, mass weighted overdensities have an increased 
dynamical range, in particular at the highest densities. As we will 
see in Sec. 13.51 this procedure, although physically motivated, 
can introduce some spurious signal, mainly induced by the mass 
of the galaxy around which the overdensity is computed. 

Another estimate of the high density environments in 
which galaxies reside c an be obtain e d by selecting optical 
groups, as described in iKnobel et al.1 (l2009h. o r X -ray ones 
dFinoguenov et al.l l2007h Finoguenov et al. l2009h ; low density 
environments can be tracked by isolated galaxies defined using 
their Voronoi volumes, as in llovino et a l. (2009). The two de- 
terminations of the environment are in a fairly goo d agreement, 
consi dering the differences of the involved scales (Kova cTet al.l 
120091) . 

In the following we use as reference the 5th nearest neigh- 
bour estimator (hereafter 5NN) of the density field, which repre- 
sents a good compromise between the smallest accessible scales 
and the reliability of the overdensity values. In this approach, 
tracer galaxies are considered within ilOOOkms -1 from the cen- 
tral galaxy and counted, after distance sorting, until their num- 
ber becomes larger than 5, considering also the fractional con- 
tribution from objects with photometric redshifts, derived from 
their likelihood function. Overdensities are then computed at the 
position of each galaxy in the spectroscopic sample, consider- 
ing also the contribution to the number or mass density of the 
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galaxy itself. We checked that the same qualitative trends of the 
GSMFs analysed in the following are present also when consid- 
ering other estimators. 



within 0.6" (~ 0.5%). The final sample contains 8450 galaxies 
with redshifts between 0.01 and 2 and 7936 in the redshift range 
where the following analysis is carried out, z - 0.1 - 1. 



2.5. Galaxy type classification 

Galaxy types can be classified in a multitude of ways, using their 
rest-frame colours, their SEDs, their spectroscopic features, their 
structural parameters and their morphologies, all of them deriv- 
able with different methods. Different classifications map dif- 
ferent physical properties. For instance, the rest-frame colour 
U - B and the galaxy SED are used as a proxy of the star for- 
mation activity and history, the morphology is an indicator of 
the dynamical state, and the two are not completely overlapping 
(iMignoli et alJl2009h . 

Even if COSMOS offers a profusion of methods to group 
galaxies, we have chosen to use only two types of classification: 
photometric and morphological. 

The photometric type is defined by SED fitting over the 
optical m agnitudes, a s sumin g as reference the same templates 
used by lllbert et all d2006l) : the 4 locally obs erved CWW 
dColeman et al.l 19801) and two starburst SEDs from lKinnev et alj 
dl996l) . extrapolated at UV and mid-IR wavelengths. These six 
templates are then interpolated to obtain 62 SEDs and opti- 
mized with VVDS spectroscopic data. The SED fitting, a y 1 
minimization performed with the code ALF (lllbert et al.l 12005; 
IZucca et all 120061 |2009l) . gives in output the best fit. Galaxies 
are then classified in 4 types, closely corresponding to colours 
of ellipticals (type 1), early spirals (type 2), late spirals (type 3) 
and irregular and starburst galaxies (type 4). We group photo- 
metric types 2, 3 and 4 in a single class, to explore in a simple 
way the evolution of the early and late-type bimodality. 

As m orphological c l assific ation we adopted the one pre- 
sented in IScarlata et al.l d2007l) : the availability of deep F814- 
band ACS@HST imag es over the whole COSMOS field 
dKoekemoer et alj 12007) allows a good determination of the 
structural parameters on which the morphology derived with 
the software ZEST (Zurich Estimator of Structural Types, 
IScarlata et alj|2007h is based. The method is a PCA analysis us- 
ing estimates of asymmetry, concentration, Gini coefficient, M20 
(the second order moment of the 20% brightest pixels) and el- 
lipticity. The morphological classes are the following: early type 
(type 1), disk (type 2, with an associated sub-classification rang- 
ing from to 3 representing the "b ulgeness", derived from the n 
Sersic indices, Sar gent et alj [2007) and irregular galaxies (type 
3). Adopting the same line of reasoning used for the photomet- 
ric types, we grouped morphologically classified galaxies in two 
broad classes, with early type including classes 1 and 2.0, i.e.,, 
ellipticals and bulge dominated galaxies. 

3. Mass Functions 

3.1. The sample 

Not all the spectroscopic redshifts have the same level of reli- 
ability, as explained in Sec. I2.ll The sample we used includes 
only the galaxies with flags corresponding to most secure red- 
shifts, i.e., starting from flag = 1 in case of agreement with pho- 
tometric redshifts. In detail, we excluded from our sample broad 
line AGNs (~ 1.8% of the statistical sample), stars (~ 5.9%), 
objects with less than 5 detected magnitudes available to com- 
pute the SED fitting (~ 1.7%) and objects for which the ground 
photometry can be affected by blending of more sources, as de- 
rived from the number of ACS sources brighter than / = 22.5 



3.2. Statistical weights 

To take into account that the observed galaxies are only a frac- 
tion of the total number of possible available targets with the 
same propertie s we need to apply statistical weight s to each ob- 
served object dZucca et al.lll994t Hibert et al]|2005l) . To this aim 
we computed the weight w, for each galaxy in our sample as the 
product of two factors connected to the Target Sampling Rate 
(TSR) and to the Spectroscopic Success Rate (SSR). Here we 
outline the basic prin ciples on wh i ch the computation is based, 
referring the reader to lZucca et al.l d2009l) for further details. 

The TSR is the fraction of sources observed in the spectro- 
scopic survey compared to the total number of objects in the 
parent photometric catalogue from which they are randomly ex- 
tracted. In t he case of zCOSM OS the VMMPS tool for mask 
preparation dBottini et al.ll2005b has been set in such a way that 
the objects have been randomly selected without any bias. A dif- 
ferent treatment has been granted to compulsory targets, i.e., ob- 
jects with forced slit positioning: they have a much higher TSR 
(~ 87%) compared with the "random" sample (~ 36%). The as- 
sociated weight is wJ SR = 1 /TSR. 

The SSR represents the fraction of observed sources with a 
successfully measured redshift: it is comprehensibly a function 
of apparent magnitude, being linked to the signal-to-noise ratio 
of the spectrum, and it ranges from 97.5% to 82% considering 
the brightest and faintest galaxies, respectively. The weight de- 
riving from the SSR is wf SR = 1 /SSR. 

The SSR is not only a function of magnitude, but also of 
redshift, since the spectral features on which the redshift mea- 
sure relies can enter o r go out the observed wavelength win- 
dow dLillvet al.ll2007b . Therefore, the redshift distribution of 
the measured redshifts can be different from the real one; it is 
possible to take into account of our lack of knowledge on the 
fail ed measure s by u sing photometric redshifts. Hence we used 
the lllbert "e t al. (2009a) release of z p h t and computed the SSR in 
Az = 0.2 redshift bins. Moreover, the characteristic emission or 
a bsorption lines are different for different galaxy types, as shown 
in lLillv et all (120091) . We further split the computation of SSR in 
each redshift bin separating red and blue galaxies, selected on 
the basis of their rest-frame U — V colour. The so-called sec- 
ondary targets, i.e., objects in the parent catalogue, but fallen 
in the slit by chance, have been considered separately: they are 
characterised by a lower SSR because they are often located at 
the spectrum edge or observed only at their outskirts. We com- 
puted and assigned the final weights w,- = wJ SR x w^ SR consider- 
ing all the described dependencies. 

3.3. Mass Function methods 

To compute the GSM Fs we adopt th e usual non-parametric 
method 1/V max dAvni & Bahcalllf jgSO). from which we derive 
the best fit Schechter function (ISchechterl 1 1 9761) . The observ- 
ability limits inside each redshift bin, z m j n and z max , have been 
computed for each galaxy from its best fit SED. 

As in iPozzetti et al.l d2009l) . we estimated the parametric fit 
of the GSMFs with both a single Schechter function, like in most 
published results, and with the sum of 2 Schechter functions, 
that seems to be statistically required by the data at least in the 
lowest redshift bins. We adopted the formalism introduced by 
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Baldrylitj4]d2004, 2006) using a single M" to limit the number 
of free parameters: 



M_ 

M* 



(1) 



Up to now the need of modeling a faint-en d upturn has been 
mainly pointed out in LF studies, both in field (IZucca et al, 119971) 
and in rich clusters (Popess oet alj 120071 iJenkins et alj 12007). 
The departure of the GSMF from a single Schechter function 
at sm all st ellar masses has be en noticed by iBaldry et al.1 ({2006. 
120081) and lPanteret~ai](l2004 on SDSS. At higher redshifts an 
a posteriori look at the published GSMFs often reveals such an 
upturn. 

We called A4 m i n the lowest mass at which the GSMF can 
be consid ered reliable and not affected by inc ompleteness on 
M/L (see lllbert et al.l2004tlPozzetti et al.l2007l). A complete de - 
scription of the procedure can be found in lPozzetti et al. (2009). 
Basically, the aim is to recover the stellar mass up to which all 
the galaxy types significantly contributing to the GSMF can be 
observed. We derived such a value in small redshift slices con- 
sidering the 20% faintest galaxies, i.e., the ones contributing to 
the low-mass end of the GSMF. For each galaxy of this subsam- 
ple we computed the "limiting mass", that is the stellar mass that 
the object would have had at the limiting magnitude of the sur- 
vey, log A1ii m = log M + 0.4(1 - 22.5). For each redshift bin we 
define as minimum mass the value corresponding to 95% of the 
distribution of limiting masses and we smooth the Almin vs z re- 
lation by means of an interpolation with a parabolic curve. We 
adopt as minimum stellar mass up to which we can reliably com- 
pute the GSMF in each considered redshift bin the Almin at the 
lowest extreme of the interval, since the 1 /V max method corrects 
the residual volume incompleteness. When dealing with GSMFs 
divided per galaxy types, the minimum masses are obtained us- 
ing the considered subsample. 

It should be noticed that this limit substantially decreases 
the number of objects considered in each redshift bin to de- 
rive the GSMF. The redshift intervals [0.10,0.35], [0.35,0.50], 
[0.50,0.70] and [0.70,1.00] were chosen to contain a similar 
number of galaxies. 

3.4. The choice of the environment definition 



As menti oned in Sec. 1X41 the density field of the COSMOS 
field (see Kov ac et ail 120 09) has been reconstructed for differ- 
ent choices of filters (fixed comoving aperture or adaptive with 
a fixed number of neighbours), tracers (flux limited or volume 
limited subsamples) and weights (stellar mass, luminosity or no 
weight, i.e., considering only the number of galaxies). 

We tested the options that allow an unbiased comparison 
over the whole redshift range, from z — 0.1 to 1.0. In particular, 
we explored the 5NN estimator and the 5NN mass weighted one 
(hereafter 5NNM), both of them computed using volume limited 
tracers, with two choices of luminosity limits: Mb < -20.5 - z 
(bright tracers) and Mb < -19.3 - z (faint tracers), where Mb is 
the absolute magnitude in the B band computed with ZEBRA. 
The absolute magnitude cut has been derived considering the 
distribution of absol ute magnitudes v s redshift, the so-called 
Spaenhauer diagram ( Spaenhauer 1978), and the evolution of the 
parameter M* B of the LFs dZucca et alj|2009l ). The requirement 
of two different limits arises from the rareness of bright tracers 
at low redshift and the incompleteness of faint tracers at high 



redshift; for this reason the two overdensity estimates cannot be 
computed over the whole redshift range, but only at [0. 1 , 0.7] 
and z - [0.4, 1 .0] for faint and bright tracers respectively. 

3.4.1 . The effect of environment tracers on GSMF 

Two issues are troubling the study on the evolution of GSMFs 
as a function of environment and must be solved: 1. we have 
to understand if the 5NNM estimator is a better tracer of the 
environment, as presumed from the theoretical point of view; 2. 
we have to be sure that the use of two different tracers, e.g., with 
a change at z = 0.7, does not introduce a spurious signal that can 
be misinterpreted as an evolutionary trend. 

To answer both questions we used as test case the redshift 
interval [0.4,0.7], where all the estimates are available, and we 
computed the quartiles of the 1+5 distribution in this redshift 
bin considering only the objects with masses larger than the min- 
imum mass. Henceforth, we will refer to the lowest and highest 
quartiles of 1 + 6 as Dl and D4 respectively. In the rest of the 
paper we will focus our study on these two extremes. 

In Fig.Q] panel (a), we show the comparison of the GSMFs, 
using a single Schechter function fit, for 5NN and 5NNM over- 
density estimators, both using the faint volume limited tracers. 
The separation of GSMFs between Dl and D4 environments is 
more prominent when considering the mass weighted estima- 
tor, because of the larger dynamical range of probed 1+6 val- 
ues. In particular, the main difference is in the massive part of 
Dl GSMF: massive galaxies in low density environments using 
5NN are moved to intermediate densities for 5NNM estimator 
because of their large stellar masses. This decreases the num- 
ber (and therefore the normalisation of the GSMF) of massive 
galaxies in low density environment when the 5NNM estimator 
is adopted. To test if this enhancement of the difference between 
the Dl and D4 GSMFs is a reliable effect, we performed the fol- 
lowing test: we removed the mass-density relation by shuffling 
the original catalogue and computing overdensities considering 
objects with their original coordinates, but assigning to each one 
the observed properties (magnitudes, stellar mass, weight) of 
the 25th following object after redshift sorting. Both 5NN and 
5NNM overdensities and their quartiles have been recomputed, 
since the shuffling implies also a change of tracers. The choice 
of the 25 object jump is a compromise between the requirements 
of preserving a similar probability to be observed at the chosen 
redshift (i.e., avoiding unphysical galaxy properties if a big jump 
in redshift is allowed) and selecting objects possibly not in the 
same structure, where we know galaxies share similar properties. 
We expected that, randomizing the contribution of stellar mass, 
the 5NN and 5NNM GSMFs would become consistent with each 
other. The comparison between GSMFs with 5NN and 5NNM 
"shuffled" overdensities is shown in Fig. [1] panel (b): a differ- 
ence between GSMFs in different environments is still present 
at high masses, now more evident in the high density environ- 
ment. 

The last test we carried out to understand this residual signal 
was the removal of the central galaxy when computing 1+6 from 
the original catalogue: the comparison of the resulting GSMFs is 
in Fig.[T]panel (c), showing now fully consistent GSMFs at high 
and low densities as defined from 5NN and 5NNM estimators. 

These tests seem to indicate that the mass weighting scheme 
assigns too much importance to the stellar mass at scales of the 
order of the galaxy itself. Then, to avoid any possible bias due 
to stellar mass over-weighting, despite its physically motivated 
link with the halo mass, we discarded the 5NNM estimator and 
performed our analysis using number weighted overdensities. 
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Fig. 1. (a) Comparison of GSMFs for environment estimates 
from 5NN and 5NNM volume limited with faint tracers: Black: 
Dl (underdense); Gray: D4 (overdense). Solid line and empty 
dots: 5NN. Dashed line and empty triangles: 5NNM. The verti- 
cal dashed line represent the value of Al m i n at z = 0.4. (b) Like 
panel (a), but 5NN and 5NNM overdensities have been estimated 
after a random shuffling of galaxy properties to remove the mass- 
density relation, (c) Like panel (a), but 5NN and 5NNM over- 
densities have been estimated without considering the properties 
of the central galaxy, (d) GSMFs for bright (M B < -20.5 - z, 
dashed lines and empty triangles) and faint (Mb < -19.3 - z, 
solid lines and empty dots) tracers using 5NN overdensities in 
the Dl (black) and D4 (gray) environments. 



Approaching the second problem, we would like to test if 
the change of the tracers at z — 0.7 may introduce some change 
in the GSMF, at risk of possible confusion with the evolution. 
We already know that the scales probed at the same 1 + 6 are 
more or less twice as large for bright compared to faint tracers 
dKovac et al.l l2009). therefore it is not possible to use the same 
1+6 threshold using faint and bright tracers. To overcome this 
problem we determined the quartiles of 1 + 6 separately for each 
redshift bin. The result of this test is shown in panel (d) of Fig. [TJ 
In the z = 0.4 - 0.7 bin, where both the tracers are available, 
the GSMFs obtained with the two tracers, with independently 
computed quartiles, are completely consistent with each other in 
under and overdense environments D 1 and D4, and therefore we 
assume we can safely compare the results at redshifts z < 0.7 
computed with faint tracers to those computed at z > 0.7 with 
the bright ones. 

3.4.2. Definition of overdensity quartiles 

As already mentioned, to trace the effect of extreme environ- 
ments on the evolution of galaxies we considered the quartiles 
Dl and D4 of the 1+6 distribution, using 5NN volume limited 
overdensities. The quartiles have been computed at each red- 
shift bin considering only the population of galaxies more mas- 
sive than Atmin of the highest redshift bin, to ensure that such 
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Fig. 2. Definition of quartiles for the 5NN estimator using vol- 
ume limited tracers: gray points represent the full sample, black 
squares the galaxies with masses above the Al m i n computed in 
the last redshift bin, horizontal segments show the values of the 
quartiles of 1 +6 computed from the distribution of these massive 
galaxies, the dashed ones indicate the median. 



a definition is not affected by the variation as a function of red- 
shift of the observable mass range, populated by different mix of 
galaxy types. The quartiles definition used throughout this paper 
is shown in Fig. [2] The median scales probed by the 5th near- 
est neighbour range from 0.87 Mpc/i^i in the D4 environment 



at low redshift to 7.57 Mpc /i^ in the Dl quartile at the highest 
redshift bin, where we are forced to use bright tracers. 

The trend toward higher values of overdensity at lower red- 
shifts is in some measure expected from the growth of structures, 
which is amplifying the dynamic range of overdensities, but this 
increase cannot be quantified through the linear approximation, 
not valid on the scales probed by our density estimates. The dif- 
ferent values of the 1+6 quartiles in the different redshift bins 
correspond to very similar scales when the same tracers are used. 



3.5. Mass Functions in different environments 

The GSMFs in the two extreme environments are shown in 
Fig. [3] the bimodal ity already visible in the global GSMFs 
( Poz zetti et al.l 120091 see also the points and lines in Fig. [3), 
with an upturn of the low-mass end around At ~ 10 9 5 Al©, is 
much exacerbated in the high density regions, at least in the two 
lowest redshift bins. We used the double Schechter function fit 
only up to z ~ 0.5, where the dip of the GSMFs falls at stel- 
lar masses larger than Al m i n - Because of the choice of the en- 
vironment definition, the normalisation of Dl and D4 GSMFs 
does not have a clear physical meaning, since the volumes oc- 
cupied by each galaxy are referred to the total volume of the 
survey and the number of galaxies in each environment is not 
1 /4 of the total sample. To have a more meaningful definition 
of the normalisation we should compute the volume occupied 
by the structures having the considered ranges in 1 +6; here 
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Table 1. Parameters of the GSMF in the low and high-density 
environments. 
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log At /Mo 
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0.35 
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+0.14 


10.53 


1.61 
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-1.12 
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0.01 
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Fig. 3. The MFs in the extreme quartiles Dl and D4 of the 5NN 
volume limited overdensities. Black: total GSMF, with 1 / V max 
dots and their Poissonian error bars and Schechter function fit 
(double Schechter function in the first two redshift ranges and a 
single one at higher redshifts). Blue: lowest 1+6 quartile. Red: 
highest density quartile. 



we are interested in the comparison of the GSMF shapes and 
hence we defer a more in-depth study of the normalisation to a 
future analysis. A striking difference in GSMFs shapes is evi- 
dent, with massive galaxies preferentially residing in high den- 
sity environments, characterised on average by a larger AT ', and 
with a steeper slope than Dl GSMFs at z > 0.35. The different 
shapes and the strong bim odality in the D4 G SMF can be ex- 
plained like the global one dPozzetti et al]|2009l) by the different 
contribution of different galaxy types, as we will see in the next 
Section. The parameters of the Schechter fits to the GSMFs are 
listed in Table Q] 



4. Evolution of galaxy types in different 
environments 

The need of a double Schechter function to fit the global and 
environment-selected GSMFs at least up to z ~ 0.5 can be linked 
to the contribution of different galaxy populations. Galaxies with 
the same luminosity can be characterised by very different Ai/L, 



Table 2. Parameters of the GSMF for the two photometric types 
in the low and high-density environments. The ratio (column still 
to be filled) is computed assuming a single Schechter function 
also for the GSMF not splitted in different classes. 







z 


a-i 


log M* /Mo 


D1T1 


0.10 


-0.35 


-0.33 


10.60 




0.35 


-0.50 


-0.17 


10.72 




0.50 


-0.70 


-0.90 


10.93 




0.70 


- 1.00 


-0.70 


10.65 


D1T2 


0.10 


-0.35 


-1.41 


10.71 




0.35 


-0.50 


-1.51 


10.81 




0.50 


-0.70 


-1.45 


10.70 




0.70 


- 1.00 


-0.84 


10.45 


D4T1 


0.10 


-0.35 


-0.03 


10.68 




0.35 


-0.50 


-0.23 


10.82 




0.50 


-0.70 


-0.28 


10.87 




0.70 


- 1.00 


-0.15 


10.94 


D4T2 


0.10 


-0.35 


-1.39 


10.92 




0.35 
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-1.43 


11.02 




0.50 


-0.70 


-1.30 


10.74 




0.70 


- 1.00 


-0.03 


10.30 



and this is why it is harder to see such a bimodal shape in LFs, 
despite this effect has been first detected in LFs, because of the 
longer history of this kind of statistical description of galaxy 
samples. 

To study the contribution of galaxies with different photo- 
metric types and morphologies in the extreme environments, we 
computed the GSMFs in Dl and D4, defined as in Sec. 13. 51 split- 
ting each sub-sample in galaxy classes. The values of A1 m i n have 
been computed separately for early/ellipticals/bulge dominated 
and late/spirals/disc dominated galaxies. These values are signif- 
icantly different, especially at low redshift, confirming the very 
different distributions of AI/L. 

The results for the contribution of different photometric 
types to Dl and D4 GSMFs are presented in Fig. [4] with the 
best fit parameters of single Schechter function fits in Table [2] 
Dividing the sample in the two broad morphological classes re- 
sults in qualitatively similar GSMFs. 

Looking at the plots in Fig. |4] it is now manifest that the 
stronger bimodality in the first two redshift bins in the D4 GSMF 
is primarily due to the larger contribution of early type galaxies. 
Just like the global GSMF, in both the considered environments 
early type galaxies are dominant at high masses (log Al/A1 ^ 
10.7), while their contribution rapidly decreases at intermediate 
masses. On the other hand, late type galaxies, being described 
by much steeper GSMFs, start to dominate at intermediate and 
low masses (log M/M Q ~ 10). 

In addition to the relative contribution of different galaxy 
types is different in Dl and D4, it is sensible asking the ques- 
tion if the shape of the GSMFs of galaxies of the same type is 
the same in different environments, i.e., if there exists a "univer- 
sal" mass function of early /late type galaxies. In Fig. [5] we show 
the comparison of early type and late type GSMFs in the two 
environments, in each redshift bin renormalised with the num- 
ber density computed for masses > Al m i n . Some small differ- 
ence on the shape of the GSMFs is present, with a trend toward 
a slightly higher density of massive galaxies in overdense re- 
gions, but the similarity of the GSMFs in all the redshift bins 
and of both the photometric types is remarkable and somewhat 
unexpected. From the plots in Fig.|5]we may argue that a quasi- 
universal GSMF characterise the different galaxy types, disre- 
garding their environment. 
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Fig. 4. Left: quartile Dl (low density environment). Right: quartile D4 (high density). Gray: total GSMF. Black: MF relative to the 
considered quartile. Red triangles and dotted lines: early photometric type galaxies. Blue squares and dashed l ines: late photometric 
type galaxies. At high masses the upper limit points show the 2cr confidence limits for detections following lGehrelsl (Il986l) . 
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Fig. 5. Left: GSMFs of photometrically early type in Dl and D4 environments, renormalised to number density = 1 for stellar 
masses > yVtmin. Right: the same for photometrically late type galaxies. Blue dotted lines and circles represent the GSMFs in 
underdense regions, Dl. Red dashed lines and squares illustrate D4 GSMFs. 



To examine the differential contribution of various galaxy 
types in different environments we can compute the evolution 
of the ratio between the GSMF of a given galaxy class and the 
global GSMFs in each environment. In Fig. [6] we show the ratio 
of 1/Vmax estimates of early type GSMF in over and underdense 
regions for the two extreme redshift bins. The trend for late type 
galaxies is reversed with respect to the one shown in the figure. 
The error bars have been computed with a Montecarlo simula- 



tion considering Gaussian distribution of errors with rms from 
Poissonian error bars from 1/V m ax method. The 16 and 84% of 
the 100,000 iterations of ratio distribution are reported in the 
plot as error bars. The vertical dashed line shows the value of 
A^min for early type galaxies in the redshift bin z = 0.7 - 1 .0. 
Despite the large error bars, Fig. [6] illustrates that at the largest 
redshifts we can probe the fractional contributions of early pho- 
tometric types to the GSMF in different environments are more 
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or less the same for Dl and D4 at all the masses we can safely 
study. On the other hand, the fractional contribution is signif- 
icantly different at low redshift, mainly at intermediate stel- 
lar masses (logyVt/yVlo 5 10.5). This trend seems to testify 
a faster growth in high density environments of the fractional 
contribution of early type galaxies. At intermediate masses the 
differences between the two extreme environments are larger: 
high stellar masses (logM/Ai Q ^ 10.7) are populated only 
by passive red galaxies in both environments, at lower masses 
(log M/Aio 55 10, where it is possible to probe them) the pop- 
ulation of late type/star forming galaxies dominates in all the 
environments. 

In a scenario consistent with these data, showing an increase 
of early type galaxies with cosmic time, intermediate mass blue 
galaxies are being transformed in more massive red galaxies, 
after quenching their star formation in a more efficient way in 
overdense regions. A possible way to quantify this different evo- 
lutionary speed is through the analysis of the evolution with red- 
shift of Across, representing the mass above which the GSMF 
is dominated by early-type galaxies. We show such a quantity 
computed from 1/V ma x points in Fig.[7]for different photometric 
types. We can see that since z ~ 1, at which the Across val- 
ues in low and high density environments were similar, the sub- 
sequent evolution produces a significant difference between the 
two Across values. The ratio of A1 cross in the highest and lowest 
redshift bins implies an evolution of a factor ~ 2 in low density 
and ~ 4.5 in high density regions. If we consider Across being 
the same in the highest redshift bin, the difference in evolution 
due to the influence of environment, computed as the ratio of the 
values A1 cross in the lowest redshift bin, is a factor ~ 3. From 
a different point of view, the plot in Fig. [7] is showing that the 
environment starts to act on the evolution of galaxies at z ~ 1, 
producing in the lowest redshift bin a delay of ~ 2 Gyr in under- 
dense regions to attain the same mix of galaxy types observed in 
high density regions. 
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Fig. 6. Evolution of the fractional contribution of the early pho- 
tometric type to the global MFs (the late type fractional contribu- 
tion is complementary to the one shown in this plot) in the two 
extreme environments. Blue lines and circles refer to low den- 
sity environment Dl (displaced by 0.02 in the abscissa to avoid 
overlapping), red lines and squares to the high density sample 
D4. Dotted lines and empty symbols represent the highest red- 
shift bin z = [0.7, 1.0], solid lines and filled points the lowest 
one, z = [0.1, 0.35]. The vertical dashed line marks A1 m i n in the 
high redshift bin (the value at low redshift is outside the plot). 
Error bars have been computed as 16 - 84% of the distribution 
of Montecarlo simulations. 



5. Discussion 

5.1. The mechanism and timescale of galaxy transformation 

Figures|6]and|7]suggest some clues on the timescale and mecha- 
nism responsible of galaxy quenching in different environments. 
We found that the evolution in the high density regions is faster 
than in low density ones, i.e., the rate of transformation in early 
photometric types is accelerated from z = 1 to low redshifts 
in overdense regions. The mechanism responsible for quench- 
ing the star formation, and then for the transformation of blue 
galaxies into passive ones, must be environment dependent. The 
physical processes operating on galaxies and transforming their 
colours and/or morphologies are manifold and can be internally 
or externally driven and gravitationally or hydrodynamically 
induc ed (see for reviews Bosel li & Gavazzil (2006; Treu et all 
120031) . The choice of possible mechanisms should exclude the 
ones acting mainly on galaxy clusters, since only a small fraction 
of our galaxies are probably located in rich clusters. We can ex- 
clude processes like ram-pressure, consisting in the gas stripping 
of a galaxy moving through a dense inter-galactic medium, with 
the effect of truncating the star formation, and harassment, i.e. 
a gravitational interaction in high velocity encounters of galax- 
ies, causing morphological transformation and bursts of star for- 
mation. Given the requirements on galaxy velocities and inter- 
galactic medium density, these processes cannot have a signifi- 
cant impact on the present study. 



Viable mechanisms in the field are instead merging and star- 
vation. Major merging processes can trigger AGN and quench 
the star formation: the fraction of pairs, related to the rate of 
merging, can be dependent on environment. The effect of merg- 
ing is to deplete the low mass end and populate the high mass one 
with elliptical morphologies. To e xplain the evolution of the den- 
sity of massive elliptical galaxies, Ilber t et aH(l200 9b) found that 
the rate of wet mergers should steeply decline at z < 1 . Limits on 
the contribution of major merging as primar y mechanism can be 
drawn from the evolution of pair fraction dde Ravel et al.l l2009b, 
who found that 20% of the stellar mass in present day galax- 
ies with log M/Mq > 9.5 has been accret ed by major merging 
events since z ~ 1) and from the GSMF dPozzettiet all [2009. 
who derived an upper limit to the fraction of major merging per 
Gyr of the order of 20 - 40% for the global population, derived 
from the GSMF evolved according to the mass growth due to 
star formation). 

Also strangulation (called as well starvation or suffocation), 
consisting of halo-gas stripping, can play a role: when the dif- 
fuse warm and hot gas reservoir in the galaxy corona is stripped 
because of gravitational interaction with low-mass group-size 
haloes or with cluster haloes at large distances from the core, 
the gas cannot be accreted anymore and the galaxy will exhaust 
the remaining cold gas through the star formation, in a timescale 
which can be instantaneous or s low, i.e., in few G yr, depending 
also on the mass of the galaxy dWolf et al.ll2.Q09l) . The result is 
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the suppression of the star formation, not followed by a morpho- 
logical transformation, explaining the possible presence of red 
spirals, even if the fading of the disc can lead to an earlier-type 
morphological classification. 



To further constrain the quest for possible transformation 
mechanisms, we computed GSMFs also div iding the sample fol - 
lowing the morphological classification by Scarl ata et al]d2007h . 
as defined in Sec. 12.51 In Fig. [8] we show the values of Across 
in the 4 considered redshift bins. This plot appears different 
from the analogous plot obtained dividing galaxies in photo- 
metric types: the values of Across are higher and their evolu- 
tion seems insensitive to the environment up to z ~ 0.4. The 
larger values of Across for morphological classification suggest 
that the dynamical transformation into elliptical galaxies fol- 
lows the quenching of the star formation. It is possible that the 
transforma tion of morphology occurs on longer timescales than 
colou r one dCapak et al.|2 007a; S mith et al.l2005l:lBamford et al.l 



120091: IWolfet alJl2009i) . as inferred also from the study of 



post-starburst galaxies selected in the same zCOSMOS sam- 
ple (Vergani et al.[2b 09) or through different evolutionary paths 
(Skibb aet al.l2008h .A more comprehensive study should be per- 
formed to investigate this point, since the larger number of pho- 
tometric early types compared with the morphological ones can 
also be due to a large fraction of dust-reddened spiral galaxies. 

In COSM OS, using photometric reds hifts between z = 0.1 
and z = 1.1, IScoville et all d2007al) and ICapak et al.l d2007al) 
find that, in the core of the large scale structures they identify, 
galaxies are characterised by earlier types, as defined by the SED 
and by the morphology, than in low density environments at the 
same redshifts, demonstrating a corre lation between galax y evo- 
lution and environment. In particular, Capa ketall (l2007al) using 
~ 33, 000 galaxies, find that the rate of transformation of spiral 
and irregular galaxies to elliptical and lenticular ones is faster 
in denser regions, whereas the star formation - density relation 
is evolving at all densities and redshifts; thus the mechanisms 
of the transformation of morphology and star formation must be 
different. 

To evaluate the uncertainties connected with this compari- 
son on photometric and morphological types, we changed the 
threshold separating elliptical galaxies from late morphologi- 
cal types: we divided the morphological class 2.1, which should 
still represent bulge dominated galaxies, following the observed 
B - z: the evolutiona ry track of the B - z colour of a galaxy Sab 
(IColeman et al.llT980T) provides a criterion to separate quiescent 
and star forming galaxies in good agreement w ith the spectral 
classification, as shown in lMignoli et al.l (12009). With this sep- 
aration, the values of morphological A4 cr0 ss become consistent 
with the photometric ones, both in the absolute values and in the 
trend with redshift. 

Both the mechanisms, gas stripping and interactions, should 
be in place to explain the suppression of the star formation and 
the morphological transformation. Those processes act in differ- 
ent timescales and are characterised by different efficiencies as a 
function of galaxy mass and environment, but it is still difficult 
to draw firm conclusions, because of the caveat associated with 
the galaxy classification. 

5.2. Comparison with mock catalogues 

We used 12 mocks COSMOS mock lightcones 
( Kitz bichler & Whit e 2007) based on the Millennium N-body 
simulation dSpringel et al.l 1200 5!). The galaxy population of 
lightco nes has then been assigned by means of semi -analytical 
recipes dCroton et al.ll200fllDeLu cia & Blaizot 2007). The final 
catalogues are the same described in lKnobel et a l. (2009), who 
used them to test the group finder algorithm. 

In this section we used 5NN flux limited 1+6 estimate of 
the environment and the rest-frame colour B-Ito separate early 
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Fig. 10. Left: quartile Dl (low density environment). Right: quartile D4 (high density). Points refer to the observed quantities, lines 
to the GSMFs derived from the mock catalogues. Black points and solid lines: GSMFs relative to the considered density quartile, 
renormalised to the same integral at log Ml M & > 10.5. Red triangles and dotted lines: galaxies with B - I > 1.15. Blue squares 
and dashed lines: galaxies with B — I < 1.15. 



and late type galaxies, to be able to compare the same quanti- 
ties in observations and mocks. Even though at the lowest stellar 
masses the mock catalogues can suffer colour incompleteness, 
this does not affect our analysis, since we limit our compari- 
son to the higher masses probed in the zCOSMOS. In Fig. [9] 
we compare the high and low-density GSMFs in the observed 
sample and in the 12 averaged mock catalogues; to avoid nor- 
malis ation uncertainties due to cosmic variance dMeneux et alj 
120091) we decided to renormalise GSMFs, in such a way that 
observed and mock GSMFs of the overdense regions have the 
same integral at masses larger than 10 10 5 Al in all the redshift 
bins. The most evident feature of the observed GSMFs, namely 
the bimodality of the GSMFs in overdense regions at low red- 
shift, is not reproduced by semi-analytical models. To explore 
what causes the failure of the semi-analytical models (SAMs) 
in reproducing observations, we separated red and blue galaxies 
adopting the threshold B - 1 = 1.15, that corresponds to the loca- 
tion of the dip of the colour bimodality, obtaining the GSMFs in 
Fig. [10] Considering first the low density environments, SAMs 
produce too many blue galaxies at intermediate and especially 
at high masses in all the redshift bins, and consequently also a 
too low density of red galaxies, in particular at 10 10 - 10 11 A1 . 
This can be ascribed to an inefficient suppression of the star for- 
mation in absence of external drivers, as in the case of sparse 
environments. Instead, in the high density regions the most vis- 
ible difference is the excess of production of low and interme- 
diate mass red galaxies (< 10 10 Ale) in SAMs with respect to 
the observed ones in the lowest redshit bin, where the probed 
mass range is wider. This last comparison reflect the problem of 
overquenching of satellites in the SAMs we used, producing too 
many small red galaxies: a too efficient strangulation produces 
an instantaneous shut down of the star format i on when a galaxy 
enters in a halo (se e Weinmann et alj 120061: iFont et alj [2 008; 
iKang & van den Boschll2008t iKimm et al.ll2009h iFontanot et al j 



2009, for a description of the issue and some attempts to solve 
it). 

6. Conclusions 

We computed GSMFs in different environments and studied the 
relative contribution of different galaxy types and their evolu- 
tion. Our main results are: 

1. the b imodality seen in the global GSMF dPozzetti et alj 
2009) up to z ~ 0.5 is considerably enhanced in high den- 
sity environments; a sum of two Schechter functions is thus 
required to reproduce the observed non-parametric estimates 
of the GSMF; 

2. the bimodality is due to the different relative contribution of 
early and late-type galaxies in different environments, each 
contribution being reasonably well represented by a single 
Schechter function; 

3. the shapes of the GSMFs of the different galaxy types and 
their evolution are quasi-universal, i.e., only the normalisa- 
tion can be significantly different in the extreme environ- 
ments we considered; 

4. the evolution of the fractional contribution of different 
galaxy types to the environmental GSMF appears to be a 
function of the overdensity in which the galaxies live, con- 
sistent with an accelerated trend of downsizing in overdense 
regions; 

5. the evolution of the crossover mass for late and early photo- 
metric type galaxies suggests a faster transition rate in over- 
dense regions, with galaxies in low-density regions experi- 
encing the same evolutionary path as the analogous galaxies 
in overdense environments with a delay of ~ 2 Gyr accumu- 
lated from z ~ 1 to z ~ 0.2; 

6. the environment starts to play a significant role in the story 
of galaxies at z S= 1 ; 
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7. the timescales for quenching of star formation and for mor- 
phological metamorphosis are different in different environ- 
ments; tentatively, the crossover mass considering morpho- 
logical classification suggest that the morphological transfor- 
mation is a delayed process compared to the colour change; 

8. SAMs fail in different ways as a function of the environ- 
ment: GSMFs computed from mock catalogues show an un- 
derestimate of the number of red massive galaxies in low 
density environments, probably due to an inefficient internal 
mechanism suppressing the star formation at relatively high 
masses; in high density regimes the overquenching problem 
of satellites in SAMs causes an excess of red galaxies at in- 
termediate and low masses. 

As a consequence of the remarkable difference in the shape 
of the GSMFs in under and overdense regions, we can infer that 
all the galaxy properties depending on mass will also depend on 
environment by virtue of the GSMF environmental dependence, 
as shown in the case of the colour-density and morp hology - 
density relations (Cucciati et al. 2009; Tasca et al.ll2009l) and of 
the AGN fraction dsilverman et al.ll2009l) . 

The nature vs nurture debate is an inextricable problem, 
since the mass of a galaxy, often referred to as its nature, is a 
strong function of the environment. A more relevant issue is the 
understanding of the mechanisms producing the observed evo- 
lution of galaxies and their transition from late to early type in 
different environments. 

Future investigations w ill also concern the imp act of merging 
in differ ent environments dde Ravel et al.l 12009a; Kampczyk et 
al.|2009) and of the role of the dark matter halo mass function in 



different environments (e.g. lAbbas & S heth 2007) in the galaxy 
formation efficiency. 
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